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Recently, a more than 2a discrepancy has been observed between the well measured inclusive 
value of V u b and the predicted value of V u b from the unitarity triangle fit using the world 
average value of sin 2(3. We attempt to resolve this tension in the context of grand unified 
SO (10) and SU(5) models where the neutrino mixing matrix is responsible for flavor changing 
neutral current at the weak scale and the models with non-proportional A-terms (can be realized 
simply in the context of intersecting D-brane models) and investigate the interplay between the 
constraints arising from B s d-B s ,d mixings, 6k, Br(r — > ^7), Br(/i — > e^) and a fit of this new 
discrepancy. We also show that the ongoing measurement of the phase of B s mixing will be 
able to identify the grand unified model. The measurement of Br(r — > erj) will also be able to 
test these scenarios, especially the models with non-proportional A-terms. 



1 Introduction 



Recent measurement of B s -B s oscillation p.J not only can examine the Kobayashi-Maskawa 
theory |2j, but also can probe the existence of new physics such as super symmetric (SUSY) 
models. Accurate measurements of the mass differences for Bd-Bd arid B s -B s , time dependent 
CP asymmetry of B^ — > J/ipK decay [3J, and Cabibbo mixing can determine the unitarity 
triangle. Since one angle and two sides are determined from these measurements, one can 
predict the remaining angles and the side. The experimental measurements of the remaining 
two angles still have large errors, while the measurement of the remaining side (naively |V^&/Vd,|) 
has become accurate recently, especially from inclusive B decay data. One can thus compare 
the experimental measurement of \V U b\ with the value from the global fit of the CKM parameters 
01 El E] and it turns out that there exists a more than 2a discrepancy between the experimental 
measurement and the unitarity prediction of IK^/V^I [U This discrepancy may be an 
indication of new physics. 

The SUSY models are the most attractive candidates for new physics. The gauge hierarchy 
problem can be solved and a natural aspect of the theory can be developed from the weak scale 
to the ultra high energy scale. In fact, the gauge coupling constants in the standard model 
(SM) can unify at a high scale using the renormalization group equations (RGEs) involving the 
particle contents of the minimal SUSY standard model (MSSM), which indicates the existence 
of grand unified theories (GUTs). While the experimental result for gauge unification is success- 
fully explained in the gauge sector, the flavor sector remains a problem in SUSY theories. The 
SUSY extension of the SM needs SUSY breaking terms, and there can be a lot of parameters 
in the flavor sector in general. The SUSY breaking scalar masses in general can induce large 
flavor changing neutral currents (FCNCs), which contradict experimental results. In order to 
avoid the SUSY FCNCs, the masses of SUSY scalar particles (squarks and sleptons) need to 
be degenerate especially for the first two generations as long as the SUSY particles are lighter 
than 2-3 TeV [ZJ. 

The flavor degeneracy of SUSY breaking scalar masses is often assumed at a cutoff scale 
such as the GUT scale or the Planck scale. The degeneracy can be realized when the Kahler 
metric is universal in flavor indices. However, even if the degeneracy is realized at a scale, the 
RGE evolution still can induce flavor violation for squarks and sleptons at low energy. The 
flavor violation induced by RGEs is suppressed by loop factors, and thus, can satisfy the current 
FCNC experimental bounds. The ongoing and future experiments may reveal the existence of 
small flavor violation, and therefore the footprint of GUT or Planck physics may be obtained 
since the pattern of flavor violation depends on the unification of matters, contents of heavy 
particles, and the nature of SUSY breaking. 
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In the scenarios where a small amount of flavor violation is generated at low energy through 
RGEs, the seeds of flavor violation are implanted in the Yukawa couplings and SUSY breaking 
scalar trilinear couplings (A-terms) . The up- and down- type Yukawa matrices are not diagonal- 
ized simultaneously, and thus, flavor mixings are generated in charged current interaction. The 
FCNCs via RGEs can originate from the mixing matrices characterized by the CKM (Cabibbo- 
Kobayashi-Maskawa) and the MNSP (Maki-Nakagawa-Sakata-Pontecorvo) matrices for quarks 
and leptons, respectively. In the MSSM, the induced FCNCs in the quark sector are not large 
since the CKM mixings are small. On the other hand, sizable FCNC effects can be generated in 
the lepton sector since two neutrino mixings in MNSP matrix are large due to the observation 
of solar and atmospheric neutrino oscillations. Thus, a testable amount of flavor violating lep- 
ton decay can be observed [8]. In grand unified models, as a consequence of the quark- lepton 
unification, the large neutrino mixings can also generate flavor violation in the quark sector [9]. 

There could be another source of flavor violation. In the minimal supergravity mediated 
SUSY breaking, the A-term couplings are proportional to the Yukawa couplings, and therefore 
the effects are already discussed above. In general, the A-term couplings can possess non- 
proportional terms which can be new sources of flavor violation. This can happen even if the 
Kahler metric is flavor universal. Actually, the A-term includes the following three parts: 1) 
universal term which originates from local SUSY breaking, 2) non-universal term which comes 
from non-trivial Kahler connection (namely, non-canonical Kahler metric), 3) non-proportional 
term which is generated when the Yukawa coupling is a function of moduli fields which acquire 
non-zero F-components. The A-term coupling in first and second parts can be proportional 
to the Yukawa coupling when the Kahler metric is flavor universal to make the SUSY break- 
ing scalar mass degenerate. In the third part, the A-term couplings are proportional to the 
derivative of the Yukawa couplings in moduli fields, and the derivative of the Yukawa matrix 
is not necessarily proportional to the Yukawa coupling itself. Thus, it is possible that A-terms 
contain the seeds of flavor violation which does not depend on the CKM and MNSP mixings, 
even if the SUSY breaking masses are degenerate at the string scale. 

In this paper, we investigate the prediction of flavor violation when the unitarity relation 
of V u b and CP asymmetry of Bd — > J/tpK is modified by the SUSY FCNC contribution. We 
assume that the SUSY breaking scalar masses are degenerate at the unification scale due to 
the nature of Kahler metric of the matter sector. The sources of flavor violation induced via 
RGEs can be the Yukawa couplings and the A-terms. We will consider the following two typical 
cases when the FCNCs are induced in SUSY breaking scalar mass matrices via RGEs. In the 
first case, the FCNCs originate from large neutrino mixings, and in the second case, the scalar 
trilinear couplings are not proportional to Yukawa couplings. 
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In the first case, we consider the SU(5) and SO(10) grand unified models. In a SU(5) 
model, the right-handed down-type quarks and the left-handed lepton doublet are unified in 
one multiplet. Thus, a sizable flavor violation can be generated in the SUSY breaking mass 
matrix for the right-handed down-type squarks. In a SO(10) model, all the matter fields are 
unified in one multiplet, and thus, both left- and right-handed squark mass matrices can have 
sizable flavor mixing. We will see that SO(10) models can have larger FCNC effects in B-B 
mixings compared to SU(5) models. 

In the second case, we will investigate the intersecting D-brane model [TUl [TT] where the 
non-proportional A-terms have been realized explicitly. In this scenario, the SUSY breaking 
scalar masses are degenerate [12] . In a simple intersecting D-brane model, the Yukawa matrices 
are rank 1 matrices plus small corrections p3J[ll|, and consequently, the non-proportional part 
of the A-term has a simple structure [15J. 

In both cases, the important constraints come from flavor violating lepton decays, fi — > 
e7 and r — > /ry, mass difference of B s -B s , and CP violation in K-K mixing. We will see 
that the experimental constraints restrict the parameters in neutrino mixings and the SUSY 
particle spectrum. Once the discrepancy between the unitarity prediction and experimental 
measurements is explained by the SUSY contribution, the phase of B s -B s mixings [T6J [17] 
and flavor violating lepton decays in these models get constrained. We will investigate the 
predictions of the models. It is interesting to note that the phase of B s -B s mixing is being 
measured by the time dependent CP asymmetry of B s — > J/ipifi decay and the CP asymmetry 
of semileptonic B s decay, and we already have a value for this phase from D0 [18] . 

This paper is organized as follows. In section 2, we will go through the unitarity relation of 
V u b and CP asymmetry of B& — > J/ ipK to understand the 2a discrepancy between the unitarity 
prediction and experimental measurement of \V u b\. In section 3, we will study the origin of 
flavor violation from the neutrino mixings in GUT models, and investigate the constraint and 
the consequences of the modification of unitarity relation by SUSY contributions. In section 4, 
we will investigate the non-proportional A-term in the context of intersecting D-brane models, 
and see the prediction of the models when it explains the 2a. Section 5 is devoted to conclusions 
and discussions. 

2 Unitarity prediction of \V u b\ 

The recent accurate measurement of B s -B s mass difference, as well as B^-B^ can determine 
one side of the unitarity triangle. In addition, one angle obtained from the CP asymmetry 
of Bd — > J/ipK and one side of the unitarity triangle determined by the Cabibbo mixing 
are accurate quantities. Therefore, the unitarity triangle can be determined accurately. The 
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remaining side and the angles are predicted if there is no effect from new physics. 

The unitarity of the CKM matrix gives rise to the following equation. The unitarity condi- 
tion, 

v ud v: b + v cd v; b + v td v* b = o, (i) 

can be rewritten to obtain the following relation, 



\v cd v; b \ 2 + \v td v; b \ 2 



2\V cd V; b V td V t t\cos(3 : 



where (3 
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The recent measurement of B s -B s oscillation indicates |Vtd/Vt s | = 0.206to;oog pp. Using world 
average of sin 2(3 obtained from B d — > J/ipK 



w 



cd | 



0.2258, and \V 



cb\ 



sin 2(3 = 0.674 ±0.026, 
(41.6 ± 0.6) x 10~ 3 [4], we obtain 
\V ub \ = (3.49 ±0.17) x 10" 3 (unitarity). 



(4) 



(5) 



The \V ub \ obtained above can be compared to the combined data analysis \ V ub \ = (3.5 ±0.18) x 
10~ 3 obtained by UTfit and CKMfitter |3 El I2D 122] ■ We note that the \V ub \ prediction from 
the unitarity triangle is insensitive to | V^^/ V* s | error since a = 2 = eng{V ud V t * d /V td V t * b ) is close 
to 90°. The prediction for a from the unitarity triangle using two sides \V cd \, \V td /Vt s \ and an 
angle (3 is a = (93.2 ± 5.7)°. The sin 2(3- V ub relation is plotted in Figure 1. 
The experimental measurement of \V ub \ [H [20] is 



\V ub \ 
\V ub \ 



(4.49 ± 0.19 ± 0.27) x 10^ 3 (inclusive), 
(3.84; 



-0.67\ 
-0.49 ) 



x 10 (exclusive). 



(6) 
(7) 



The tension between the experimental data and the unitarity prediction of \V ub \ is significant 
especially for the inclusive B decay data which has less theoretical error compared to the 
exclusive case. The unitarity condition does not agree with the inclusive data within 99% 
confidence level. The data from the exclusive decay B — > -niv includes large uncertainty from 
lattice calculation. The statistical error of the inclusive data is reduced in recent experiments 
at Babar and Belle [20, 23J, and the inclusive data has become more reliable. 

The discrepancy between the unitarity condition and the experimental measurement may 
indicate new physics. We consider a scenario in which the unitarity condition, Eq. ([!}), is really 
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Figure 1: The la range of the unitarity relation between \V u b\ and sin 2/3 (solid lines). The 
central value of experimental measurements is plotted at x. We also plot the experimental 
region for 68%, 90%, and 99% confidence levels assuming that the errors are Gaussian. 

true, but the SUSY contribution modifies the sin 2/3 measurement by Bd — > J/ipK. We will 
study different types of boundary conditions (SU(5), SO(10) GUTs or non-proportional A- 
terms at the unification scale) that can explain the discrepancy without contradicting other 
experimental results. 

3 Origin of FCNC in GUT Models 

In SUSY theories, the SUSY breaking terms can be the sources of flavor violations. In general, it 
is easy to include sources of flavor violation by hand and the discrepancy between the unitarity 
condition and the experimental measurement can be fitted since the SUSY breaking masses 
with flavor indices are parameters in the model. However, if these parameters are completely 
general, too much FCNCs are induced. Therefore, as a minimal assumption of the SUSY 
breaking, universality is often considered, which means that all the SUSY breaking scalar 
masses are universal to be Wq, and the scalar trilinear couplings are proportional to Yukawa 
couplings (the coefficient is universal to be Aq) at unification scale. In this case, the angles in 
the unitarity triangle are not changed. Since the measurements of the V u b are obtained from 
tree-level processes, the SUSY particles do not contribute to its determination. Therefore, we 
need to include flavor violating sources at the unification scale to explain the discrepancy. 

In GUT models, the flavor violating sources in SUSY breaking parameters can be induced 
from neutrino Dirac Yukawa couplings, Y U 5N C H 5 [HIS]- Since the left-handed lepton doublet, 
L, and the right-handed down-type quarks, D c , are unified in 5 multiplet, the SUSY breaking 
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mass matrix for D c is corrected by colored-Higgs and right-handed neutrino loops. As a result, 
the flavor violation in the quark sector can be generated from the neutrino coupling. The 
contribution is naively proportional to Y V YJ. Our purpose in this section is to search for a 
solution of the V u b-sm2/3 discrepancy using the neutrino mixings in the context of a GUT 
scenario. 

We will work in a basis where the right-handed neutrino Majorana mass matrix, Mr, and 
charged-lepton Yukawa matrix, Y e , are diagonal, 

Mr = diag (Mi, M 2 , M 3 ). (8) 

The neutrino Dirac Yukawa coupling matrix is written as 

Y v = V[Y^V^\ (9) 

where V£ R are diagonalizing unitary matrices. We note that V£ corresponds to the (conjugate 
of) MNSP neutrino mixing matrix, £7mnsp ; i n type I seesaw, m l l ght = Y^M^Yj (H®) 2 [24], up 
to a diagonal phase matrix when V£ = 1 (identity matrix), which we will assume for simplicity. 
Through RGE, the off-diagonal elements of the SUSY breaking mass matrix for the left-handed 
lepton doublet gets the following correction 

* -^( 3 ™o + <) E(nm*ln^, (10) 

k k 

where M* is a cutoff scale and the SUSY breaking parameters are universal. Neglecting the 
threshold of the GUT and the Majorana mass scales, we can write down the boundary conditions 
as 

( I h 

1-kV£\ k 2 \V?\, (IV, 



Mi = M% = M\ = m\ 



where k ~ (^ diag )| 3 (3 + A 2 / m 2 ) /8rr 2 In M,/M GVT , and k 2 ~ Am 2 JAm 2 atm M 2 /M 3 . We pa- 
rameterize the unitary matrix V£ as 



/ i(ai—S) \ / „e „e „e „e „e e iS 

I ( \ I c 12 c 13 A '12 L 13 5 13 e 

i 12 c 23 c 12 i 23 i 13 e c 12 c 23 6 12 i, 13 6 23 e c 13 i 23 

*12*23 c 12 c 23 i '13 e L 12*23 i, 12 ± '13 c 23 e c 13 c 23 



V 



(12) 

where sf- and cf- are sin and cos of mixing angles In the limit fci^ — * 0, «i and «2 are the 
phase of the 13 and the 23 element of M|. Since we are assuming that V£ — 1, #12 and #23 
correspond to solar and atmospheric neutrino mixings, respectively, which are large. On the 
other hand, sf 3 is bounded by CHOOZ experiments, s^ 3 <^ 0.2 [23]. The SUSY breaking mass 
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for 10 multiplet (Q, U c , E c ) is also corrected by the ( colored- )Higgsino loop, but it arises from 
CKM mixings and the effect is small. So, we assume that the boundary condition at the GUT 
scale for 10 multiplet is 

Ml Q = M| = M? c = M| c = ml 1, (13) 

neglecting the small contribution arising from the CKM mixings. The boundary conditions, 
Eqs.f fTTfTBl . are typical assumptions in the case of SU(5) GUT. The Yukawa coupling matrices 
for up- and down-type quarks and charged-leptons are given as 

Y u = V qeL V^ KM Y^P u V uR , (14) 
Y d = V qeL Yt g P d V q T eR , (15) 
Y e = Y e diag P e , (16) 

where Y u are real (positive) diagonal matrices and P u ,d,e are diagonal phase matrices. In a 
minimal SU(5) GUT, in which only H 5 and H$ couple to matter fields, we have V u r = VckmV^ l , 
VqeL = V qe R = 1, and Y d diag = Y" e diag . We may consider non-minimal SU(5), but we assume 
V u r, V qe L, VqeR — 1- Otherwise, the simple relations between the flavor violations in the quark 
and lepton sectors which we see below will be lost. 

All matter fields are unified in the spinor representation 16 in SO(10) models. Since the 
right-handed neutrino is also unified to other matter fields, the neutrino Dirac Yukawa coupling 
does not have large mixings (i.e. V£ ~ 1) in a simple fit of the Yukawa couplings. In this 
case, the proper neutrino masses with large mixings can be generated by the type II seesaw 
mechanism [26J, in which the interaction term \fLLA L (where A L is a SU(2) L triplet) induces 
the light neutrino masses, m^ ght = f(A L ). Due to the unification under SO(10), the left-handed 
Majorana coupling, /, is unified to the other matter fields, and the off-diagonal terms in the 
sparticle masses are induced by loop effect which are proportional to ff>. Neglecting the GUT 
scale threshold, we can write the boundary condition in SO(10) as 



1-k 16 u( 1 k 2 jE/tj, (17) 
where k 16 ~ 15/4 (/ 3 d 3 ag ) 2 (3 + Al / ml) /8n 2 In M J M Gl]T , and k 2 ~ Am 2 sol /Aml tm in this case. 



M 2 W = ml 



Note that the parameters K\§, k\ )2 are of course different from those given in Eq. pTT) using 
the set-up for type I seesaw, but we use the same notation to make the description simple. 
The unitary matrix U is the (conjugate of) MNSP neutrino mixing matrix up to a diagonal 
phase matrix, which is parameterized in the same way as Eq. ffT2l) . The Yukawa couplings are 
also given as Eq. fll4lll5|fT6l) . If we do not employ 120 Higgs fields, the Yukawa matrices are 
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symmetric, and thus, V u r = VckuY^ l , V qe R = V qe L- The unitary matrix V qe L is expected to be 
close to 1 if there is no huge fine-tuning in the fermion mass fits. 

Now let us consider the necessary conditions to solve the discrepancy between the unitarity 
condition and experimental measurements using the boundary conditions, Eqs. (1111)131) in SU(5) 
and Eq. ffTTl) in SO(10) respectively. Due to the smallness of k\$ and s\ 3 , the 12 and 13 elements 
of M| and Mf 6 are expected to be smaller than the 23 element. In the limit of ki^, sf 3 — > 0, 
only the 23 element is non-zero and then sin 2/3 can not be modified. In this limit, the B s -B s 
mass difference Am s is modified, and therefore |V^/Vt s | prediction is changed, and one of the 
sides of the unitarity triangle is modified. However, if |V^/V^| is modified to obtain larger 
\Vub\ from the unitarity relation, CP violation in K-K mixing, ex, will not fit well with the 
experimental range. Thus, we need to modify sin 2/3 itself employing the finite contributions of 
the 13 element in SUSY breaking mass matrix. If at least one of k 2 and sf 3 is non-zero, the 13 
element is generated. Then the 12 element of SUSY breaking mass matrix is also generated in 
general. However, if the 12 element is generated, due to the quark-lepton unification, /i — > 
process is enhanced through the chargino diagram. In order to satisfy the experimental bound 
for the branching ratio of fi — > e'-f process, Br(/i — > ej) < 1.2 x 10~ n [27] . the 12 element 
needs to be small to suppress the chargino contribution. It can happen by a cancellation when 
k\ <C &2, s\ 3 ~ &2 sin 26*12/2 and 5 ~ n. 

In order to obtain a sizable effect in sin 2/3, k needs to be large enough. However, if k is 
too large, r — > process may exceed the experimental bound, Br(r — > yy) < 4.5 x 10~ 8 [28J 
for large #23 mixing. A large value of tuq ~ 1 TeV can avoid the excess of r —>■ /ry due to the 
following reason. The off-diagonal elements of SUSY breaking mass matrices are proportional 
to Mq, whereas the diagonal element of squark matrix at the weak scale is insensitive to m 
when mo ^ 500 GeV for gluino mass ~ 1 TeV, since RGE contribution from the gluino is 
large. Consequently, the flavor violations from SUSY contributions in the quark sector become 
maximal for mo ~ 1 TeV. On the other hand, the diagonal elements of slepton mass matrix are 
sensitive to m , and thus, lepton flavor violating processes are suppressed for large m . 

Now let us study different results that can be produced by type I seesaw for SU(5) boundary 
condition, Eq. ( fTTyTBl . and type II seesaw for SO(10) boundary conditions, Eq.tflTl). If the 
boundary condition is flavor universal, which means k = 0, the chargino contribution is the 
dominant SUSY contribution to the B ds -B ds mixing amplitude, M 12 (-Bd jS ). Using a general 
parameter space for the soft SUSY breaking terms, the gluino box diagram dominates the SUSY 
contribution. The gluino contribution (divided by the SM contribution) can be written naively 
in the following mass insertion form 

ttM* * a [(S d LL )l + (8UI] - b (S d LL U8 d RR ) 3i , (18) 



S 



where i — \ for mixing, i = 2 for B s mixing, a and b depend on squark and gluino masses, 
and Sf L RR = (Mj} LL}RR /m 2 (m is an averaged squark mass). The matrix Mj is a down-type 
squark mass matrix (Q, D cJ( )Mj(Q' f , D C ) T at weak scale in the basis where down-type quark 
mass matrix is real (positive) diagonal. When the squark and the gluino masses are less than 
1 TeV, a ~ 0(1) and b ~ 0(100). We also have contributions from 8f R , but we neglect them 
since they are suppressed by (mb/msusy) 2 - 

Using the SU(5) boundary condition, Eq. ( TTTffTBl . we find that 5 RR can be large due to large 
neutrino mixings, but 5f L can not be very large since it is induced from CKM related mixing. 
On the other hand, using the SO(10) boundary conditions we find that both <5^ L and 5 RR can 
be large. As a result, the SO(10) boundary condition can generate larger SUSY contribution 
in B-B mixings rather than the SU(5) boundary condition for the same parameters for k, k\ i2 
and mixing angles 0^ since b 3> a in Eq. ljTB]) . 

The B-B mixing amplitude is given as 

M 12 (5) fu11 = M 12 (B) SM + M 12 (5) SUSY . (19) 

The angle 2(5 for unitarity triangle is the argument of Mi 2 (Bd) SM in the particle data group 
phase notation of CKM matrix. The observed angle 2(3 eS by Bd — > J / ipK decay is the argument 
of M 12 (5 ( i) fu11 . The SUSY contribution of B^-B^ mixing is naively proportional to (s^ 3 ) 2 . As 
was noted, it is necessary that s\ 3 ~ fc 2 sin#i 2 /2 to suppress fi — > e7 process. In the case 
of type I seesaw in SU(5) model, k 2 is a free parameter, and thus, s{ 3 is free. On the other 
hand, using type II seesaw in the SO(10) model, k 2 is almost determined by the ratio of mass 
squared differences for solar and atmospheric neutrino oscillations. As a result, the suppression 
of ji — > e^f process indicates s\ 3 ~ 0.02. Due to the smallness of s^ 3 , the modification of sin 2(3 
will be small for type II seesaw SO(10) boundary condition. For type I seesaw SU(5) boundary 
condition, on the other hand, the modification can be larger than the SO (10) case since s^ 3 can 
be sizable ~ 0.1, though the modification is smaller if we use the mixing angle ~ 0.02. 

The absolute value of Mi 2 (-E>d(s)) SUSY is almost determined by k, s^ 3 (s| 3 ), and squark, 
gluino masses, whereas, the argument of M 12 (S) SUSY can be anything since ai. 2 in Eq. (1121) and 
phases P^ jU in Yukawa matrices in Eq. fll4lll5p are all free. The minimal value of sin2/? eff for a 
given size of SUSY contribution can be obtained by choosing the phases. The branching ratios 
of lepton flavor violating processes do not depend on the phases since the branching ratios are 
proportional to the squared absolute values of the amplitudes of SUSY contributions. 

We plot the relation between the minimal value of sin 2/3 eS and the branching ratio of r — > /i7 
decay in Figure 2. In the plot, we use mo = 1.2 TeV, = 300 GeV for gaugino mass at 
GUT scale, and tan/3# = 10 for a ratio of vacuum expectation values of up- and down-type 
Higgs fields. In the SU(5) case, one can tune k 2 for a given s^ 3 to suppress Br(/i — > ej) to 
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be less than 10 -11 . With type II seesaw in the SO(IO) model, k 2 is fixed by observation of 
neutrino oscillation and the suppression of Br(/x — > ej) implies sf 3 ~ 0.02. Since there are 
also neutralino contributions for both left- and right-handed muon decay, one can not cancel 
Br(/i — > ej) completely. In the range where sin 2/3 is in agreement with experiment, we obtain 
Br(/i — > e^y) > 10~ 12 — 10~ 13 . In the plot, we show the cases s{ 3 = 0.08, and 0.15 for SU(5), 
and sf 3 = 0.022 for SO(10). As was expected, one can read off that the modification of sin2/? efr 
can become large for larger s{ 3 for the same value of Br(r — > /ry). 

So far we have obtained the parameter region which satisfies the observed sin 2/? eff and the 
inclusive data of \V u b\, as well as the bounds on the leptonic flavor violating decay in a grand 
unified scenario. We will next study the consequence of the fit in the phase of B s -B s mixing. 
Although the measurements of mass differences, AM^ S = 2|Mi 2 (.Bd iS ) full |, are in good agreement 
with SM prediction, there can still be room for SUSY contribution due to the following two 
reasons. One reason is that there exists large uncertainties ~ 20% in the lattice calculation 
for the decay constants fs da and the bag parameters Bs da - The SM prediction of the ratio 
AM S / AMd can be more accurate since the ratio for fs and Bb have less uncertainty. The 
other reason is that we still have an ambiguity for the phase of B s -B s mixing. For instance, 
one can satisfy |Mi 2 (£? s ) full | = |Mi 2 (-B s ) SM | by choosing a phase of Mi 2 (B s ) SUSY as long as 
|Mi 2 ( J B s ) SUSY | < 2|Mi 2 (.B s ) SM |. As we have noted, the phase of Mi 2 ( J B s ) SUSY is completely 
free for the boundary conditions in GUT models. The SM prediction of the phase of B s -B s 
oscillation (2/5, = -arg M 12 {B S ) SM ) is about 2° (sin2/5 s ~ 0.04). If the SUSY contribution 
modifies sin2/3 eff and resolves the tension of the sin2f3-V u b unitarity relation, the phase of 
B s -B s (2f3f = -arg M 12 (B s Y un ) must be different from the SM prediction. 
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Figure 3: Random dot plot for sin2/? eff and 4>b s for \V u b\ = 0.0041 using the SO(10) boundary 
condition. The shaded area shows the la range of sin 2/3 for the world average [20J . 



It is convenient to define the deviation from SM prediction by two real parameters Cb s and 
4> Bs as 

c 2i4> Bs = M 12(^) fUU (2Q) 

UBs€ ~ M 12 (5 s )sm ' ^ 

If there is no SUSY contribution, Cb s = 1 and 4>b s =0. In this notation, /3| ff = j3 s — (pB s - We are 
interested in the value of 4>b s in the range where sin2/3 e fj is in agreement with the experimental 
data for inclusive V^. In the case of SU(5) boundary condition, the value is 4>b s ~ ±3°, namely 
sin2/?| ff ~ +0.14, —0.07. On the other hand, in the case of SO(10) boundary condition, 4>s a can 
be sizable. We show <pB s and sin2/5 e g in SO(10) type II seesaw by random dot plot in Figure 
3. We choose the same parameters as in the plot for Figure 2. Constraint arising from the 
Br(r — > /ry) < 4.5 x 10~ 8 rules out larger values of \4>b 3 \ ^ 35°. We also filter out the points by 
the constraint from K-K mixing, and AM S /AM^. Especially, points at the left-bottom corner 
in the plot are ruled out by K-K mixing data. The reason for the allowed left-top corner is that 
the SUSY contribution becomes 6^^^^ ~ -2e| M . From the fi gure, we find the preferred value to 
be <pB s ~ +(20 — 30)°, which corresponds to sin2/5| ff ~ —(0.6 — 0.9). The recent measurement 
by D0 collaboration [18] shows a large negative central value for sin 2/3| ff = —0.71 with a sizable 
error and if this result holds then, definitely, the SO(10) model will be preferred. 

We summarize what we have accomplished in this section. We have considered the pa- 
rameter space where we can resolve the discrepancy between the unitarity prediction and the 
experimental measurement for the boundary conditions for SU(5) with type I seesaw and for 
SO(10) with type II seesaw. To obtain the solution we need large m ~ 1 TeV to suppress 
t — > /i7 decay. In order to suppress fi — > e^y, we need to adjust the 13 mixings in the Dirac 
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Yukawa couplings for type I seesaw and the Majorana couplings for type II seesaw. Conse- 
quently, the neutrino oscillation parameter is constrained as follows. The CP phase in neutrino 
oscillation <5mnsp needs to be close to ir. In type II seesaw case, the 13 neutrino mixing should 
be around 0.02, while in type I seesaw case, the 13 mixing is free. The phase of B s -B s is not 
same as that can be obtained in the SM. The deviation is small for the SU(5) case, and large 
for the SO(10) case. 

The crucial assumptions we made to obtain above results are that ~ 1 in the Dirac 
Yukawa coupling, Eq.([H]), in type I seesaw for SU(5), and the triplet part dominates the light 
neutrino mass rather than the type I seesaw part for type II seesaw for SO (10). If we do not 
assume them, V£ and U in the boundary conditions do not take part in the construction of 
the neutrino mixing matrix. Even in this case, we may need the relation, s\ 3 ~ k 2 sin 2^12/2, 
and 5 ~ -n to suppress ji — > decay, but the mixing angles and the phase are not necessarily 
related to the observed neutrino mixings. Moreover, since s| 3 is not necessarily large, r — > fi'-f 
can be suppressed by choosing a small s| 3 even for smaller <^ 500 GeV. If this is the case, 
depending on the s| 3 , the deviation of the B s -B s phase (pB 3 will not be large compared to the 
above. The possible feature of this case is that Br(r — > ej) may be comparable to Br(r — > /ry) 
since s\% needs to be around 0.1 to modify sin2/5 efr while we get Br(r — > cy) <C Br(r — > fij) in 
a usual scenario. 

4 Origin of FCNC Using Non-proportional A-term 

In the previous section, we have considered the scenario where the neutrino mixings are the 
origin of flavor violation at low energy, assuming the universality for SUSY breaking scalar 
masses and trilinear scalar couplings at unification scale. The SUSY breaking scalar mass and 
trilinear scalar couplings are given as 

Ka = m 2 3/2 +V - £ F M F*d M d„ In K a - a , (21) 

M,N 

A abc = £ F M [(K M - d M ln(K a - a K bl K c5 ))Y abc + d M Y abc ] , (22) 

M 

where K is the Kahler potential and F M is the F-term of moduli field M. The flavor in- 
visibleness of Kahler potential is necessary for the universality of SUSY breaking parameters. 
However, the invisibleness is not a sufficient condition to make the A-term flavor universal. 
Actually, when the Yukawa coupling depends on moduli and the F-term of the moduli is non- 
zero, the v4-terms are no longer proportional to Yukawa coupling, even if the Kahler potential 
is flavor invisible. In the intersecting D-brane models [11] , the Kahler potential can be flavor 
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invisible, but Yukawa couplings depend on [/-moduli (shape moduli), and they can generate 
flavor violation in the SUSY breaking scalar mass matrices at low energy |15j . In this section, 
we will consider non-proportional A-terms to resolve the sin 2/3- tension. 

In the intersecting D-brane models, the SUSY breaking scalar masses are flavor universal 
when the generations are simply replicated at the intersection of the branes. The Yukawa 
couplings are induced by the three-point open string scattering, and they are determined by 
the triangle area formed by branes. As we have noted, the couplings depends on [/-moduli. 
The Yukawa coupling matrices are given in a factorized form 

" - xf(C/ 1 )xf(C/ 2 ), 



(23) 



if the left- and right-handed matter fields are replicated on different tori. The couplings are 

1 (t). The coupling Xi is given as 



the function of Jacobi theta function § 



X\ : %2 : X3 = $ 



e + 




(t):0 







(t):0 



£ 




(t), 



(24) 



where e stands for a brane shift parameter, and the Wilson line phase is chosen to be zero 
for simplicity. The variable t is proportional to [/-moduli. Since the Yukawa matrices are 
rank 1, only third generation becomes massive after the electroweak phase transition. We, 
therefore, need to introduce corrections to Yukawa couplings, which can arise due to higher 
dimensional operators. If the correction SY is given as 5Y = diag (0, 0, e) in the basis where the 
light neutrino Majorana mass is diagonal, f/ e 3 (the 13 mixing in neutrino oscillation) is exactly 
zero. In the quark sector, the correction 5Y = diag (0, 0, e) leads to the relation m s /m h ~ V co . 
The first generation fermions are still massless with the above correction to the rank 1 Yukawa 
matrix and we need more corrections. Thus, U e ^ is small and relates to the size of the Cabibbo 
angle in the quark-lepton unification picture. Furthermore, tan9 so \ <> 1 for the solar mixing 
angle, and sin 2 2# atm ~ 1 for the atmospheric neutrino mixing can be satisfied naturally in the 
perturbative region of the Yukawa coupling. The details can be found in the Refs. [HI [15] . 

If the F-terms of [/-moduli, F u , are zero and the Kahler potential is flavor invisible, the A- 
term is proportional to the Yukawa couplings. We consider the case where the Kahler potential 
is flavor invisible whereas F u 7^ 0. The rank 1 matrix Y° is given as in Eq.(23). For simplicity, 
we consider that the matrix is symmetric: xf' R = Xi(U). The diagonalization unitary matrix 
of Y° is 



cos 9 S 
cos 6 a sin 9 S 
sin 9 a sin 9 S 



— sin 9 S 

cos 9 a COS 9 S 

sin 9 n cos 6L 





- sin 9 

cos 9 n 



) 



(25) 



where tan# s = x\jxi and tan# a = Jx\ + x^/x?,- We work in a basis where the light neutrino 
mass is diagonal assuming that the correction for the charged-lepton Yukawa matrix is SY e ~ 
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diag(0, 0, e). Then, MNSP neutrino mixing matrix is given as 



^MNSP —V*U , 



(26) 



where V e is a diagonalization matrix of U°(Y° + 5Y e )(U°) T . If the quark-lepton unification 
is realized, we obtain V e ~ Vckm- The solar and atmospheric neutrino mixing is close to 
6 SA respectively, but solar mixing is corrected by (V e )\2. The 13 neutrino mixing is given as 



The non-proportional part of A-term is proportional to djjY. The derivative of Yukawa 
coupling Y° is calculated as 



where y = Y^xf. The ij = 1,2) elements of dY° are zero, and thus, the 12 element of 
non-proportional the A-term is suppressed to be of the order of ~ X 2 dY^ + A 3 <9Y23 (where 
A ~ 0.2) in the basis where the charged-lepton Yukawa matrix is diagonal. This is a general 
property of the almost rank 1 Yukawa matrix, independent of any set-up. 

We note that the non-proportional pieces of the A-term can be considered in the non- 
minimal scenario of hidden sector model with flavor symmetries. However, in such a scenario, 
the non-proportional A-term is hierarchical and thus, large 13 off-diagonal elements are not 
expected. On the other hand, in our scenario of £7-moduli origin of the A-terms, we can realize 
large 13 elements of the A-term. 

In general, if the A-term has non-proportional pieces, the off-diagonal elements of SUSY 
breaking scalar mass matrices are generated via RGEs, e.g. dM^ c / \\iQ = 1/(Att 2 )AJA^ + • • •. 
Since our purpose is to modify sin 2/5, we need the 13 (31) element of the A-term for left 
(right) handed squark mass matrices in the basis where Yukawa matrices are diagonal. The 
RGE induced 13 element of SUSY breaking mass matrices can be larger when A33 {=A b ) 
is large. Therefore, a large coefficient of A-term, i.e. Aq, is preferable. The negative Aq 
can generate large contribution to sin 2/3 for gluino mass M 3 > since the RGE is given as 
dA b /d\iaQ = +8/ (3Tr)a 3 M 3 y b + • • •, where y b is a bottom quark Yukawa coupling. 

When both A 13 and A 2 3 elements are introduced at boundary condition, the 12 elements of 
SUSY breaking mass matrices are generated at low energy even if A12 is zero in the boundary 
condition. The 12 elements of SUSY breaking mass matrices are unwanted for K-K mixing 
and fi — > e-y. Thus, for our purpose to modify sin 2/?, we need to suppress A 2 3 (A 32 ) rather than 
A13 (A 31 ). Such a situation can be satisfied when d\jd a = 0. There is a solution of djjO a = 
when i < s < | where e is a brane shift parameter in Eq. (l2"l"|) . Since d\jQ a = 0, the atmospheric 



U e3 ~ (V e )l 2 /V2. 




/ 


\ ysm9 a du9 s yd v 9 { 



a 




(27) 
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Figure 4: The sin2/? eff and Br(r — > ej) relation employing non-proportional A-term. The 
dotted lines stand for the experimental bounds. 

mixing as a function of U is maximized for a given e. The predicted atmospheric mixing is 
| < sin 2 2# a < 1 and 9 a > 45°. 

For the numerical studies, we assume the following boundary condition for the non-proportional 
A-term: 



in the basis where Yukawa matrices are diagonal. The A-terms for up- and down- type quarks 



A np 's are unified, but we assume the unification in the context of G422 = SU(4) c x SU(2)lX 
SU(2)i? model construction of intersecting D-brane: The A np may originate from G422 invariant 
term, while the corrections to the Yukawa coupling breaks G i2 2- The A np is assumed to be 
symmetric just for simplicity. 

We plot sin2/? efr and Br(r — > ej) in Figure 4, varying 5 in Eq. (|28|) from to y^. In the 
plot, we use m = 200 GeV, m 1/2 = 300 GeV, tan/3 H = 10 and \V ub \ = 4.1 x 10~ 3 . When 
sga(yd) = — sgn(?/fe) for Yukawa couplings of down and bottom quarks, the SUSY contribution 
of sin 2/3 becomes negative and satisfy the observation from CP asymmetry of Bd — > J / ipK. 
The positive values of Aq do not modify sin 2(3 much as we have explained. The Br(r — > e^y) 
do not depend on the signature of A because RGE running of A T is not large. 

Since we have assumed that A12 and A23 are zero, ex and B s -B s mixing (AM, and </>b b ) do 
not get modified from the SM values. We do not get any enhancement in the branching ratios 
of [A — ► e7 and r — > fij, neither. When we introduce A23 elements, all these quantities can be 




(28) 



and charged-leptons are given as A u ^, e 



AoY U) d, e + A^ de . In general, there is no reason that 
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modified. The strongest constraint comes form ex when sin 2/3 is in the experimental range. 

Interestingly, in this scenario we have solutions for smaller values of m . Thus, the solution 
can be consistent with the muon g — 2 from e + e~ data [29] . The main feature of this solution 
is that Br(r — > ej) is larger than Br(r — > /X7). 

5 Conclusion 

The measurements of B^s-B^s mass differences, AM^ S , and the CP asymmetry of — > 
J/ipK, sin 2/3, determines the unitarity triangle, and generates predictions for other angles, 
CP violation in K-K mixing, as well as |Kife/Kb|- The predictions are consistent with the 
measurements within errors except for \V u b/V c b\. The recent development of \V u b\ measurement 
from the inclusive B decay indicates a more than 2a discrepancy compared to the value obtained 
from the unitarity relation. This may imply the presence of new physics. 

Since the inclusive B decay to measure \V u b\ is a tree level process, it does not include new 
physics. On the other hand, B-B mixing and sin 2/3 arise from box diagrams, and thus, they 
may include new physics contributions, since the new particles can propagate in the diagram 
and can shift the experimental measurements of AM Sj ^ and sin 2(3. 

We have studied whether the SUSY contributions can modify the unitarity relation to 
make the experimental measurements to fit well, and investigated the predictions for other 
measurements, such as the CP asymmetry in B s decays and branching ratios of flavor violating 
lepton decay modes. As a consequence, the modification of the unitarity relation of |T4&| and 
sin 2(3 shifts the tension to the phase of B s -B s mixing or Br(r — > e7)/Br(r — > yrf). 

We assumed that the SUSY breaking scalar masses are universal at a cutoff scale. The 
origin of flavor violations can be in the Yukawa interaction and scalar trilinear couplings. We 
considered two different set-ups to modify the unitarity relation. In the first set-up, the flavor 
violation originates from the large neutrino mixings in SU(5) and SO(10) grand unified models. 
In the second set-up, A-term couplings include pieces which are not proportional to Yukawa 
couplings in the context of intersecting D-brane models. 

In GUT models, responsible for quark-lepton unification, the neutrino Dirac couplings and 
the Majorana couplings can induce the off-diagonal elements of squark mass matrices via the 
GUT scale particle loops. In order to modify the unitarity relation appropriately, we need 
large contribution in 13 elements in SUSY breaking scalar mass matrices. However, a large 13 
neutrino mixing is bounded by experiments, since it leads to large effects in flavor violating 
lepton decay such as r — > 7x7 and // — > ej. When we suppress the branching ratio of the lepton 
flavor violating decays, neutrino mixing parameters and the SUSY spectrum is restricted. As 
a result, the CP violating phase in the neutrino oscillation needs to be close to tt when the 
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large solar and atmospheric mixings are responsible for the structures of neutrino Dirac Yukawa 
coupling in type I seesaw and for the Majorana coupling in type II seesaw. In type II seesaw 
case, the 13 neutrino mixing is restricted to be around 0.02. The squarks and the sleptons need 
to be heavy for the solution, and the universal scalar mass mo is about 1 TeV. As a consequence 
of resolving the V u b-sm2f3 tension by the SUSY contribution in the GUT models, the phase 
of B s -B s mixing becomes different from the SM prediction, /3f M ~ 1°. The prediction of the 
phase in the SU(5) model is not very different from the SM prediction, /3f u ( 5 ) ~ ±3°, while it 
becomes large in the SO(10) boundary condition to be /3f°( 10 ) ~ ±(20 - 30)°. In the SO(10) 
case, the contribution of ex is also large. In order to fit the ex data, the negative value of /3| ff 
is preferred. The D0 collaboration has reported the direct measurement of the phase, and the 
phase /3g can be a large negative value, /3| ff = (—23 ± 16)° [18]. The combined analysis of 
the CP asymmetry in the B s — > J/ip'P decay as well as asymmetry in semileptonic B s decay 
improves = (—16 ± 12)°. The measurements are expected to become more accurate soon 
and the prediction of the models can be tested. 

The v4-term can include a non-minimal term which is not proportional to the Yukawa 
coupling if the Yukawa coupling is the function of moduli and the moduli acquires a non- 
zero F-component which breaks SUSY. In the intersecting D-brane models, it is natural that 
the SUSY breaking scalar mass matrices are flavor universal while the v4-term can have non- 
proportional terms. In this model, when the atmospheric neutrino mixing angle is maximized 
as a function of [/-moduli, only the 13 and the 31 elements of the non-proportional A-term 
can be large. As a result, the phase of Bd-Bd mixing can be modified to make the unitarity 
relation consistent with the \V u b\ data. The feature of this model is that Br(r — > cy) is larger 
than Br(r — > fij). The current bounds of the branching ratios are Br(r — > /ry) < 4.5 x 10~ 8 
and Br(r -> ej) < 1.2 x 10~ 7 [28]. 

Finally, we comment on the sin2/9 eff from b — > s penguin modes such as Ba — > 4>K de- 
cay in our solutions to resolve the V / u f ) -sin2/3 tension. In GUT models, since the squarks are 
heavy ~ 1 TeV in the solutions, the contribution of the penguin diagrams are small. Thus, 
sin 2/3]^ d _+j/Tp K ~ sin 2j3 e ^ d ^ K . In our solutions using the non-proportional /L-terms, the b — ► s 
penguin effects can become larger depending on the model parameters since the squarks are 
light. However, when 23 and 32 elements of the non-proportional A-term are zero as we have 
assumed in Eq. (l28l . the b — > s penguin effects are not significant. 
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